Spin quantum Hall effect and plateau transitions in multilayer network models 
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We study the spin quantum Hall effect and transitions between Hall plateaus in quasi two-dimensional net- 
work models consisting of several coupled layers. Systems exhibiting the spin quantum Hall effect belong to 
class C in the symmetry classification for Anderson localisation, and for network models in this class there is an 
established mapping between the quantum problem and a classical one involving random walks. This mapping 
permits numerical studies of plateau transitions in much larger samples than for other symmetry classes, and we 
use it to examine localisation in systems consisting of n weakly coupled layers. Standard scaling ideas lead one 
to expect n distinct plateau transitions, but in the case of the unitary symmetry class this conclusion has been 
questioned. Focussing on a two-layer model, we demonstrate that there are two separate plateau transitions, 
with the same critical properties as in a single-layer model, even for very weak interlayer coupling. 

PACS numbers: 72.15.Rn 64.60.De 05.40.Fb 
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I. INTRODUCTION 

Universality classes for critical behaviour at Anderson tran- 
sitions are determined by the dimensionality and the symme- 
tries of the HamiltonianJ- The best known universality classes 
are the three Wigner-Dyson classes originally identified in the 
context of random matrix theory. Seven additional symmetry 
classes for localisation were recognised^ over a decade ago: 
they are distinguished from the Wigner-Dyson classes by hav- 
ing a special point in the energy spectrum and energy levels 
that appear in pairs either side of this point. In one of these 
additional classes, known as class C, properties of suitably 
chosen quantum lattice models can be expressed in terms of 
observables for a classical model defined on the same lattice. 
This mapping was originally discovered^ in the context of the 
spin quantum Hall effect (SQHE), where it relates a delocal- 
isation transition in two dimensions to classical percolation, 
also in two dimensions, for which many relevant aspects of 
critical behaviour are known exactly. 

Models in class C arise from the Bogoliubov de-Gennes 
Hamiltonian for quasiparticles in a gapless, disordered spin- 
singlet superconductor with broken time-reversal symmetry 
for orbital motion but negligible Zeeman splitting. The spe- 
cial energy in this case is the chemical potential (which we set 
to zero) and pairs of levels are related by particle-hole sym- 
metry, which has profound consequences for the influence of 
disorder on quasiparticle eigenstates. The quantum to classi- 
cal mapping provides a framework within which quasiparticle 
properties can be studied in great detail starting from a sim- 
plified description of a disordered superconductor. 

Here we use this approach to study for the SQHE an aspect 
of the plateau transition that has resisted detailed investigation 
in the context of the conventional integer quantum Hall effect 
(IQHE) belonging to the unitary symmetry class. Specifically, 
we study plateau transitions in models of n weakly coupled 
layers. For such systems, scaling ideas"* and the sigma model 
description^ lead one to expect n distinct transitions, separat- 
ing adjacent pairs of phases in which the Hall conductance 



differs by one quantum unit, irrespective of how weakly the 
layers are coupled. The same scaling ideas have other impor- 
tant consequences. In particular, they suggest a scenario for 
the disappearance of the IQHE as magnetic field strength is 
reduced, in which extended states responsible for plateau tran- 
sitions levitate^'^ in energy, and they are input for construction 
of the global phase diagram,** in which the Hall conductance 
of adjacent IQHE phases again differs by one quantum unit. 
Alternative types of behaviour have also been proposed, in- 
volving direct transitions between phases with Hall conduc- 
tance differing by multiple quantum units^^ and there has de- 
veloped quite an extensive literature on the subject, reviewed 
in Ref. 10. Attempts' '"'^ to distinguish between these dif- 
ferent possibilities using numerical simulations are hampered 
by the fact that in the most interesting regimes - weak mag- 
netic field or weak interlayer coupling - the localisation length 
is never short, making asymptotic behaviour hard to reach. 
Even in one of the simplest settings, involving a two-layer 
IQHE system as a representation of a spin-degenerate Landau 
level, the existence of two transitions has been inferred only 
rather indirectly. By contrast, we show in the following for 
the SQHE that the mapping to a classical description allows 
simulation of sufficiently large systems that the behaviour ex- 
pected from scaling and the sigma model can be revealed in 
considerable detail. 

The mapping^ between a single-layer network model for 
the SQHE and classical percolation generalises to all models 
in the same symmetry class that share a set of key features.— 
A variety of physical quantities of interest for localisation (al- 
though not all) can be determined in this way, including the 
two-terminal conductance of a finite sample, which will be 
our main tool. The generalisation, however, relates the quan- 
tum problem to a classical one involving interacting random 
walks rather than percolation, so that while much is known an- 
alytically about percolation, simulations are required to study 
properties of the classical walks. Since these classical simu- 
lations are much less computationally intensive than a direct 
study of the quantum problem, much larger system sizes are 



accessible. Here we exploit the classical mapping to study the 
conductance of quasi two-dimensional SQHE systems. Sim- 
ilar conclusions to the ones we present have been suggested 
in earlier work,'^ but from a much more restricted range of 
system sizes. Our results for quasi two-dimensional sys- 
tems are complementary to recent work on the metal-insulator 
transition in a three-dimensional class C network model, in 
which the classical mapping enabled a measurement of crit- 
ical exponents with a precision unprecedented for a localisa- 
tion transition."* 

While our focus is on properties of the quantum system, 
we believe that the classical walks we study deserve attention 
in their own right. In particular, it would be interesting if ar- 
guments could be found directly for the classical problem, to 
show that the n-layer system has n transitions, and that these 
are in the same universality class as classical percolation on 
the plane. 



II. MODEL 

We study models in which quasiparticles propagate along 
the directed links of a lattice and scatter between links at 
nodes. Disorder enters the models in the form of quenched 
random phase shifts associated with propagation on links. For 
class C models, the disorder-averaged (spin) conductance can 
be expressed as an average over configurations of interact- 
ing classical random walks on the same directed lattice .^^^iii 
This relation between quantum properties and classical walks 
holds on any graph in which all nodes have exactly two in- 
coming and two outgoing links. In the classical problem the 
connection at each node is a quenched random variable hav- 
ing two possible arrangements. Incoming and outgoing links 
are arranged in pairs, and a particle passing through the node 
follows the pairing with probability p, or switches with prob- 
ability 1 — p. Given a directed graph with the required co- 
ordination, any choice of classical connections at the nodes 
separates paths on the graph into a set of distinct, closed, mu- 
tually avoiding walks. Average properties of these walks are 
calculated from a sum over all node configurations, weighted 
according to their probabilities. 

Here we consider systems formed from several coupled lay- 
ers. Each layer is an L x i square sample of the L-lattice, 
shown schematically in Fig. 1. It is characterised by the 
node probability p and in isolation has a plateau transition at 
p = 1/2. We construct n-layer models by stacking n copies 
of this lattice in register, with independent disorder realisa- 
tions in each layer. For n = 2 we couple the layers using a 
second set of nodes, located at the mid-points of the links in 
each layer This set of nodes is characterised by the probabil- 
ity pi of switching layers, and the model is symmetric under 
pi — > 1 — pi. At pi =0, and also at pi — 1, the system sep- 
arates into independent layers, both with a plateau transition 
at p = 1/2; we are concerned with behaviour as a function 
of p for < pi < 1, and in particular whether two coupled 
layers exhibit two transitions at separate values of p. We also 
examine, though in less detail, a three-layer system, consider- 
ing only a parameterless form of interlayer coupling. This is 
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FIG. 1 : (Color online) The L-lattice: circles represent nodes at which 
links are connected as indicated with probability p, and in the oppo- 
site sense with probability 1 — p. Left and right panels indicate the 
boundary conditions applied to obtain the quantities we refer to re- 
spectively as the longitudinal conductance and Hall conductance. In 
both cases, current leads are attached to the left and right edges. For 
calculation of the longitudinal conductance, periodic boundary con- 
ditions are applied to the top and bottom edges to form a cylinder. 
For the Hall conductance, reflecting boundary conditions are applied 
to the top and bottom edges. 

We calculate the two-terminal conductance between oppo- 
site, open faces as the average number of classical trajecto- 
ries connecting these two faces. We use two different types 
of boundary condition in the other direction, as illustrated in 
Fig. 1. With periodic boundary conditions, so that the sample 
is a cylinder, the two-terminal conductance of a system de- 
scribed by a constant local conductivity tensor would be sim- 
ply the longitudinal conductivity. For that reason we call the 
average conductance in this geometry the longitudinal con- 
ductance. Alternatively, using reflecting boundary conditions, 
the two-terminal conductance within a Hall plateau is deter- 
mined by the number of edge states. We therefore call the 
average conductance in this geometry the Hall conductance, 
even though its value between Hall plateaus depends on both 
components of the conductivity tensor. Within the framework 
of the quantum-to-classical mapping we use, the disorder- 
averaged spin conductance G{p, L) of the quantum system 
is given (in units of h/in) by the average of the number of 
classical paths from a specified open face to the other 

Our simulations use system sizes L of between 500 and 
5000 lattice spacings. For the largest system, we average up 
to 10^ disorder reaUsations. Earlier work^^ was limited to 
L < 80. 



III. CONDUCTANCE AND SPIN QUANTUM HALL 
TRANSITIONS 

We first study the two-layer system at interlayer coupling 
Pi — 1 /2. In Fig.|2]we show the behaviour of the longitudinal 
and transverse conductances as a function of the intralayer pa- 
rameter p for two system sizes. Two transitions are apparent, 
at p ~ 0.43 and p 0.57, separating three phases charac- 
terised by quantised Hall conductances of 0, 1 and 2 units. 
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The accurate quantisation of the Hall conductance in the cen- 
tral phase is striking. In terms of classical walks, it arises 
because almost all realisations that contribute to the average 
have exactly one extended trajectory at each reflecting edge. 
Other trajectories in this phase are typically much shorter than 
sample size, since the longitudinal conductance is very small. 
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FIG. 2: (Color online) Conductance of the two-layer system at pi — 
1/2 as a function of p for system sizes L = 1000 (dashed lines) and 
4000 (continuous lines). Curves with steps (blue) represent the Hall 
conductance, obtained with reflecting boundary conditions. Curves 
with peaks (red) are the longitudinal conductance, obtained with a 
cylindrical sample. The inset shows the scaling of the longitudinal 
conductance near the first peak as a function of {p — pc)L^^^ for 
L = 3000, 4000 and 5000. 



These transitions in the two-layer system are expected to 
have the same critical behaviour as for a single layer, which 
maps to that of the percolation transition. As a first test, we 
plot in the inset to Fig. |2] the longitudinal conductance near 
the peak at p ~ 0.43 as a function of {p — Pc)L^^'^ for three 
system sizes L = 3000 (blue dots), 4000 (green triangles) and 
5000 (red squares). The good overlap of the data for differ- 
ent sizes shows that the widths of the peaks in longitudinal 
conductance scale as L^^/", where v — 4/3 is the correlation 
length exponent for percolation. The width of the steps in Hall 
conductance follows a similar behaviour 

As a second and more precise check of universality, we 
examine the peak value of the longitudinal conductance. In 
a single layer this can be deduced from the crossing proba- 
bilities for percolation clusters in an annulus. We obtain 
values for these from Ref. |2l|: the probability of a single 
crossing in a square sample with cylindrical boundary con- 
ditions is 0.357369 . . . and the probability of two crossings is 
0.002018 . . ., while higher order crossings are neglegible at 
our precision. So the exact value of the critical conductance 
in this geometry from percolation theory is 0.361404 .... To 
test our numerical approach, we first calculate the maximum 
value of the longitudinal conductance for a single layer using 
a sequence of system sizes at the exact critical probability for 
percolation, p = 1/2. We find that deviations of these maxima 
from the large system limit decrease roughly proportionally 
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FIG. 3: (Color online) Conductance of a single-layer system at crit- 
icality, p = 1/2, versus L~^. The large dot on the vertical axis 
corresponds to the exact value. 



to L^^, so we extrapolate to infinite size by plotting maxima 
as a function of L^^. We show the result in Fig. [3] A lin- 
ear fit yields the intercept 0.36135 ± 0.00010, confirming the 
reliability of the approach. We also obtain the maximum con- 
ductance for a bilayer system, but this time we do not know 
the exact critical probability for percolation so we have to cal- 
culate the conductance for several probabilities in the critical 
region and fit the results with a Gaussian to extract the maxi- 
mum value at each system size. The deviations of this maxi- 
mum value from the large system limit do not scale as either 
range of system sizes considered. To 
extract a limiting value of the critical conductance for the bi- 
layer system, we therefore first plot the difference between the 
bilayer maximum conductances and the critical conductance 
given by percolation theory, as a function of L on a double- 
logarithmic scale. The result is shown in the inset to Fig. |4l 
the data lie close to a straight line with slope x = 0.662. We 
then plot in the main part of Fig.|4]the conductance as a func- 
tion of L^^. This figure shows that the conductance extrapo- 
lates to a value close to the one from percolation theory, and 
we obtain 0.3613 ± 0.0002. The uncertainty is higher than for 
a single layer because the critical value of p is not known ex- 
actly and because the system sizes employed are necessarily 
smaller. 

We have also studied the conductance for the three-layer 
model described above. Results are shown in Fig. 5. As 
expected, three distinct transitions separate four phases, and 
adjacent phases have Hall conductances differing by a single 
unit; the critical points are at p ~ 0.37, 0.5 and 0.63. 



IV. PHASE DIAGRAM FOR A TWO LAYER SYSTEM 

The results presented in the previous section are for maxi- 
mal coupling between layers, and it is interesting to examine 
how behaviour changes as this coupling is reduced. In partic- 
ular, a key question is whether the degenerate transition occur- 
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FIG. 4: (Color online) Maximum value of the conductance of a two- 
layer system as a function of L""'®''^. The large dot on the vertical 
axis corresponds to the exact value for percolation. Inset: maximum 
conductance, with respect to the percolation value, versus L on a 
double-logarithmic scale. 
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FIG. 5: (Color online) Conductance as a function of -p for the three- 
layer system defined in the text, with system sizes L — 1000 (dashed 
lines) and 4000 (continuous lines). Curves with steps (blue) represent 
the Hall conductance, obtained with reflecting boundary conditions. 
Curves with peaks (red) are the longitudinal conductance, obtained 
with a cylindrical sample. 



ring atp — 1/2 for n uncoupled layers is split into n distinct 
transitions by any non-zero interlayer coupling, or whether a 
step in Hall conductance of more than one unit persists for 
small couplings. We focus on properties of the two-layer sys- 
tem as a function of the interlayer coupling pi, and calculate 
the longitudinal conductance rather than Hall conductance be- 
cause it is easier to locate the transition probabilities using 
this quantity. For all system sizes and values of pi at which 
there are two clear peaks in the longitudinal conductance as 
a function of p, we determine peak positions p^ by fitting the 
conductance of one of them to a Gaussian inp — pc- At each 



value of pi we make an extrapolation of these peak positions 
to infinite system size, linearly in L^". We find empirically 
that a = 2 is the best fitting exponent in all cases considered, 
but the results do not depend much on this value. We consider 
sizes L from 500 to 5000. In large systems it is possible to 
identify two distinct transitions even at very weak interlayer 
coupling. For example, withpi = 5 • 10"'' we can distinguish 
two peaks for sizes greater than 1000, and forpi = 10"'' we 
can discriminate for sizes greater than 500. The extrapolated 
values of the critical probability for the transitions are shown 
inFig.|6] 
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FIG. 6: (Color online) Phase boundaries for spin quantum Hall states 
in a two-layer system, as a function of the intralayer parameter p and 
the interlayer parameter pi . Points are from extrapolations described 
in the main text; the continuous curve is a guide to the eye. Inset: 
same data on a double logarithmic scale. 

To gain further insight, we analyse in detail the shape of the 
phase boundary at small interlayer coupling. In this regime 
Pc is close to 1/2 and if two distinct transitions persist for all 
non-zero pi, we expect a relation of the form 

PI cx Ipc - l/2r. (1) 

In the inset to Fig. 4 we show pi versus pc — 1/2, on a double 
logarithmic scale. The straight line is a linear fit to the four 
points with smallest pi. Its slope is 7 = 2.004 ± 0.002. The 
fact that the phase boundary is accurately described by Eq. ([TJ 
is good evidence for the correctness of the quantum Hall scal- 
ing flow diagram of Ref. 4 when applied to the SQHE. It 
is also evidence that in the system studied there is no direct 
transition between phases with Hall conductance differing by 
more than one unit. 

The value 7 = 2 can be understood by the following ar- 
gument. Consider first a single-layer system, and let ^ (p) be 
the coiTelation length for classical walks: the typical diame- 
ter of the largest closed loops. These loops are the hulls of 
percolation clusters, and for large ^ (p) their arc-length varies 
as ^(p)'''" with dh — 7/4. Moreover, from the mapping-' for 
a single layer to classical percolation, we know that ^(p) di- 
verges as ^(p) oc |p — 1/21"" with = 4/3 when p ap- 
proaches 1/2. Next examine the probability in a two-layer 
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system that a pair of such loops, one from each layer, are 
coupled. This probability is expected to be of order one for 
p — Pc- It is made up of the product of three factors: the 
length ^(p)'^'' of one loop, the density ^(p)'*'*^^ of the other 
loop and the probability pi that a given pair of links in equiv- 
alent positions in the two layers are coupled. We therefore 
expect pi (X ^(pc)^"^'^'' and hence pi oc |p — 1/2 p'^^'*''^^) = 
\p— 1/2 p. It is interesting to note the difference between this 
result for the SQHE and the equivalent one for a two-layer 



IQHE system, in which (taking over the notation we have 
defined for the SQHE) pi cx \pc — 1/2|''. In the language 
of the quantum localisation problem, this difference arises be- 
cause the density of states vanishes at the mobility edge for 
the SQHE but is finite for the IQHE. 
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